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Abstract 


The  purpose  of  this  research  is  to  determine  the  funda- 
mental performance  improvement  of  an  optimum  array  detector 
versus  a conventional  beamformer  detector  in  a jamming 
environment.  The  fundamental  performance  improvement  is 
based  on  total  knowledge  of  both  the  desired  signal  and  the 
jammers,  and  their  respective  locations  in  space.  The 
jammers  are  modeled  as  a colored  noise  component  in  the 
binary  hypothesis  detection  problem.  Mathematical  tract- 
ability  is  achieved  by  considering  distributed  measurements 
in  space  and  time  across  the  array.  Conditions  are  identified 
which  allow  determination  of  the  necessary  eigenvalues  and 
eigenfunctions  by  inspection  of  the  propagation  functions  of 
plane  wave  jammers. 

Performance  improvement  is  given  for  several  jamming 
scenarios:  two  mutually  orthogonal  jammers,  single  jammer 

(two  spatial  dimensions),  band  limited  jammer,  and  single 
Jammer  (one  spatial  dimension).  The  optimum  detector  performs 
respectively  68dB,  60dB,  19.5dB,  and  20dB  better  than  the 
beamformer  for  the  jamming  scenarios  addressed  and  similar 
parameters  for  each  case.  The  maximum  performance  improve- 
ment for  the  two  mutually  orthogonal  jammers  occurs  when 
twenty-five  percent  of  the  signal  power  projects  in  the 
direction  of  each  jammer.  Maximum  improvement  for  a single 
Jammer  occurs  when  fifty  percent  of  the  signal'  power  projects 
in  the  direction  of  the  Jammer.  Performance  Improvement  in 


the  presence  of  a spatially  bandlimited  jammer  is  due 
primarily  to  the  degradation  of  the  beamformer's  perfor- 
mance. For  the  scenarios  examined,  the  optimum  detector 
provides  significant  performance  improvement  over  the 
conventional  beamformer  detector. 
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PERFORMANCE  IMPROVEMENT  OF  A CLASS 
OF  ARRAY  PROCESSORS  IN  A JAMMING  ENVIRONMENT 

I.  Introduction 


Recently  much  research  has  been  devoted  to  the  subject 
of  multiple  sensor  or  array  type  antenna  systems.  This 
interest  has  been  generated  by  the  ability  of  the  adaptive 
array  antenna  systems  to  provide  rapid  inertialess  scanning 
for  high  speed  angular  coverage,  their  large  power  handling 
capability,  and  their  usability  in  environments  where  actual 
movement  of  the  antenna  is  difficult  or  impossible  (Gallop, 

1971:  2).  When  using  an  array  system  (or  any  other  antenna 
system),  the  goal  is  to  process  the  measurement  provided  by 
the  array  to  extract  information  contained  in  the  received 
signal  as  accurately  as  possible.  Heretofore,  the  design  of 
a radio  frequency  (RF)  communications  receiving  system  has 
been  divided  into  two  segments:  the  antenna  system  and  the 
communications  processor.  Typically,  each  of  these  segments 
was  designed  and  optimized  separately.  These  two  optimally 
designed  segments  were  then  united  to  form  a single  system 
which,  due  to  the  ad  hoc  nature  of  this  design  approach, 
may  have  resulted  in  an  overall  sub-optimal  receiver  processor. 

Such  a design  approach  has  been  motivated  by  several 
factors.  First,  the  disciplines  of  antenna  design  and  processor 
design  were  treated  as  separate  and  quite  distinct  areas  of 
academic  interest  and  application.  Second,  the  optimum  proces- 
sor has  sometimes  led  to  computations  which  were  unrealizable 
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from  a hardware  point  of  view.  This  limitation  has  been 
somewhat  overcome  recently  by  technology  with  the  advent  of 
practical  applications  of  digital  processing  techniques, 
surface  acoustic  wave  devices,  and  charge  coupled  devices 
(Dudgeon,  1977:98;  Melen  and  Buss,  1977:327-328). 

With  the  ability  to  design  and  potentially  build  optimum 
receiver  systems  at  hand,  the  increased  performance  gained 
from  the  optimum  system  needs  to  be  compared  with  the 
performance  of  the  conventionally  designed  system.  This 
comparison  then  becomes  a useful  tool  in  assessing  the  cost 
relationships  related  to  the  two  design  approaches. 

Depending  upon  the  type  of  signal  to  be  detected  by 
the  antenna  array  processor  system,  different  criteria 
exist  for  the  performance  comparison  of  the  optimum  and  con- 
ventional array  systems.  For  the  case  of  a deterministic 
signal,  the  comparison  of  the  output  signal  to  noise  ratios 
(SNR)  of  each  system  results  in  a useful  measurement  of 
performance  improvement  of  the  optimum  antenna  array  proces- 
sor over  the  conventionally  designed  system  (Van  Trees,  1968:99). 
This  ratio  of  SNR  is  the  method  adopted  to  evaluate  the  funda- 
mental performance  improvement  of  the  optimum  array  processor 
in  this  thesis. 

Some  work  has  been  done  in  the  area  of  determining  the 
performance  of  optimum  adaptive  array  receivers.  Gallop  and 
Nolte  addressed  the  performance  of  an  array  receiver  for  a 
signal  of  unknown  spatial  location  in  spatially  uncorrelated 
Gaussian  noise.  They  demonstrated  a trade-off  among  array 
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parameters  and  signal  parameters  relevant  to  the  performance 
of  the  detector  (Gallop  and  Nolte,  1974:429-435).  In  his 
dissertation.  Gallop  concluded  that  he  was  unable  to  calcu- 
late performance  analytically  even  from  the  closed  form 
expression  for  discrete  sampling  by  a two  or  three  sensor 
array  for  a signal  at  an  unknown  location  (Gallop,  1971:73). 
Adams  and  Nolte  studied  the  performance  of  the  array  for  a 
signal  of  known  location  but  an  uncertain  waveform  in  the 
presence  of  spatially  uncorrelated  noise.  They  found  that 
the  detection  performance  was  equal  to  the  performance  of 
the  scalar  processor  with  input  signal  to  noise  ratio 
increased  by  a scalar  factor  (Adams  and  Nolte,  1975:656-669). 
Hodgkiss  and  Nolte  addressed  the  performance  degradation  for 
the  array  processor  when  uncertainty  exists  in  the  direction 
of  the  signal,  noise  source,  or  both  (Hodgkiss  and  Nolte, 
1976:605-615).  All  of  the  performance  results  obtained  by 
the  previous  studies  included  the  effects  of  unknown  signal 
properties  (direction,  power,  etc.)  as  well  as  uncertainty  in 
the  noise  source  properties.  The  inclusion  of  these  uncertain 
properties  clouds  the  issue  of  fundamentally  how  much  better 
the  optimum  array  receiver  does  perform  over  any  other  array 
receiver  in  the  presence  of  jamming. 

Pasupathy  developed  a performance  comparison  between 
an  optimum  and  conventional  array  processor  using  a distri- 
buted (spatially  continuous)  measurement  to  simplify  the 
mathematics  associated  with  a discrete  element  array.  His 
work,  however,  addressed  only  one  spatial  dimension  and  one 
noise  model  for  a passive  sonar  signal  problem  and  did  not 
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include  the  effects  of  signal  or  noise  uncertainty  (Pasupathy, 

x978: 158-164) . 

Examination  of  the  fundamental  performance  improvement 
of  an  optimum  array  processor  in  the  presence  of  spatially 
correlated  noise  (jammers)  is  the  central  issue  of  this  thesis. 
The  approach  to  be  followed  extends  Pasupathy 's  work  to  two 
spatial  dimensions  and  examines  several  different  jamming 
environments.  The  performance  of  the  optimum  processor  is 
considered  a benchmark  since  the  signal  of  interest  is  assumed 
to  be  fully  known  and  thus  represents  the  most  improvement 
which  the  optimum  array  processor  can  provide  in  the  jamming 
environment . 

The  following  chapters  will  address  the  representation 
of  the  jamming  signals  as  spatially  colored  noise.  A review 
of  detection  theory  will  be  followed  by  the  performance 
improvement  results  of  the  optimum  array  processor.  Con- 
clusions drawn  from  these  results  and  any  recommendation  for 
additional  research  will  be  presented  in  the  final  chapter. 


* 

\ 
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II . Statistical  Representation  of  Jamming  Signals 

This  chapter  addresses  four  topics  in  developing  the 
statistical  representation  of  jamming  signals.  These  topics 
are:  (1)  scalar  field  theory,  (2)  fields  as  random  processes, 

(3)  second  moment  description  of  jammer  fields,  and  (4)  series 
representation  of  jammers. 

A.  Scalar  Field  Theory 

Before  looking  at  any  representation  of  a jamming  signal 
it  is  necessary  to  establish  a notation  convention  which  is 
based  on  scalar  field  theory.  Consider  first  the  scalar 
field  given  by: 

u(t,x,y,z)  = A(t,x,y,x)cos[2Tifot-<Kt,x,y,'z)]  (1) 

which  is  a function  of  both  time  and  space  coordinates.  This 
expression  can  more  easily  be  written  as: 

u( t ,x,y , x ) = Re{U(t,x,x,z)exp[-j2nfQt]}  (2) 

where 

U(t,x,y,x)  = A(t,x,y,z)exp[j<J>(t,x,y,z)]  (3) 

U(t,x,y,z)  is  called  the  complex  envelope  of  the  scalar  field. 
This  representation  is  appropriate  for  both  electric  and 
magnetic  fields.  For  the  purpose  at  hand,  the  representation 
is  of  an  electric  field  with  units  of  volt s/me ter/ /ohm.  It  is 
easy  to  see  that  the  complex  envelope  contains  all  the  infor- 
mation of  the  field  centered  at  fQ. 

Scalar  field  theory  and  the  use  of  the  complex  envelope 
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permit  a straightf orward  means  of  representing  a signal 
impinging  on  an  array  surface.  Fig.  1 shows  the  signal 
originating  from  a point  source,  propagating  through  space 
and  striking  the  array  surface.  When  the  distance  between 
the  array  surface  and  the  point  source  is  large  enough  then 
the  wavefront  striking  the  array  surface  is  a plane  wave 
(Gagliardi  and  Karp,  1967:13).  For  ease  of  notation  it  is 
convenient  to  consider  the  field  at  a specific  value  of  the 
z coordinate,  thus  reducing  the  complex  envelope  represen- 
tation to  two  spatial  dimensions.  With  this  notation  in 
mind,  the  complex  envelope  of  the  plane  wave  striking  the 
array  is  given  by: 

U(t,x,y)  = A(t)exp[j>(t)]exp  [j  (vxx+vyy)]  (M 

where 

2ir  Q 
v = — yCOS0 

X A 

and 


where  A is  the  wavelength  of  the  plane  wave  given  by: 

A = (5) 

where  c is  the  speed  of  light  and  f is  the  frequency  of  the 
wave.  A(t)  is  the  time  varying  amplitude  of  the  plane  wave 
and  ijj(t)  is  the  time  varying  phase  of  the  plane  wave.  When 
A(t)  and  i|>(t)  are  constant  for  all  values  of  t then  the  plane 
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wave  is  said  to  be  monochromatic.  When  the  bandwidth,  B, 
is  much  less  than  the  center  frequency,  f , the  plane  wave 
is  said  to  be  quasi-monochromatic  (Gagliardi  and  Karp, 
1976:13-17).  Otherwise  the  plane  wave  is  non-monochromatic . 
B.  Fields  as  Random  Processes 

When  the  fields  of  interest  are  unknown,  it  is  convenient 
to  model  them  as  complex  random  processes  in  time  and  space. 
The  purpose  of  this  section  is  to  develop  a notation  for 
these  random  processes. 

Consider  first  a complex  representation  of  the  field: 


U(t,r)  = UR(t,r)  +.  JUjCt,?)  (6) 

where  UR(t,r)  and  U^Ct,?)  are  the  real  and  imaginary  parts  of 
U(t,r).  The  notation  r is  the  vector  representation  of  the 
point  described  by  the  coordinates  x and  y. 

The  first  moment  of  the  random  field  is  defined  as: 


E[U(t,r)]  = E[UR(t  ,r)  ] + JEtUjU,?)] 

= MR ( t , r ) + JMjCt,?)  = M(t,r)  (7) 

Given  the  first  moment  of  the  random  field,  the  second 
moments  can  be  obtained  from  the  correlation  functions, 
E[UR(t,r)UR(t',r ' ) ] and  E[Uj  ( t ,r  )UI  ( t',r ' ) ] , and  the  cross 
correlation,  E[UR(t , r )Uj  ( t',r ' ) ] . Direct  evaluation  shows  that 
the  same  information  is  contained  in  the  complex  terms 
E[U(  t , r )U(  t’,r ' ) ] and  E[U(  t , r )U#  ( t',  r')  ] for  complex  envelopes. 

When  working  with  random  fields  typical  assumptions  are 
that  the  real  and  imaginary  parts  have  the  same  correlation 
functions  and  that  they  are  uncorrelated  and  have  zero  mean. 
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Hence , 

E[U(t,r)U*(t\  ?'  )]  = R(  t , r , t',r ' ) 

= 2E[UR(t,r)UR(t',r')]  (8) 

where  R(t,r,t',r')  is  the  correlation  function  of  the  field 
represented  by  U(t,r).  Several  cases  of  the  correlation 
function  are  of  interest  for  later  use: 

(a)  Temporally  Stationary: 

R(t,r,t',r ' ) = R(r,?,x)  (9a) 

where  t = t-t' 

(b)  Spatially  Stationary: 

R ( t , r , t’,  r ' ) = R ( r-r ' , t , t ' ) (9b) 

(c)  Coherence  Separable 

R( t , r , t ' , r ' ) = Rt(t,t')  Rg(r,r')  (9c) 

where  R.(t,t')  and  R (?,?')  are  the  temporal  and  spatial 
correlation  functions  respectively. 

C . Second  Moment  Description  of  Jammer  Fields 

The  question  of  how  to  model  a jamming  environment  can  be 
considered  from  a number  of  aspects  depending  upon  the  infor- 
mation available  about  the  jammer.  For  the  case  at  hand,  con- 
sider a set  of  point  sources  in  the  far  field  of  the  array 
of  which  one  of  the  point  sources  is  the  signal  of  interest 
and  the  other  point  sources  are  jammers  as  shown  in  Fig  2. 

Since  the  point  sources  are  in  the  far  field,  this  implies 
that  the  wave  fronts  arriving  at  the  array  are  nearly  plane 
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Figure  2.  Point  sources  in  far  field  of  array  representing  jammers  and  desired 
signal  (s). 


waves.  The  detection  problem  for  the  jamming  environment 
becomes  one  of  deciding  whether  or  not  the  desired  signal 
is  present  in  the  plane  waves  striking  the  array.  In  order 
to  solve  this  detection  problem,  it  is  necessary  to  look 
at  the  second  moment  properties  of  the  jammers. 

To  account  for  both  jammers  and  the  thermal  noise  of 
the  array  (always  present),  the  total  noise  interference  is 
considered  to  be  a spatially  non-white  Gaussian  noise  composed 
of  an  independent  white  component  which  is  both  temporally  and 
spatially  stationary  representing  the  thermal  noise  of  the 
array  and  a colored  component  representing  the  jammer  (Van  Trees, 
1968:287-288).  The  use  of  a colored  component  representation 
allows  for  the  inclusion  of  several  jammer  scenarios  without 
significantly  changing  the  detection  problem  for  the  array. 
Generally,  for  the  array  problem,  the  noise  term  is  given  by: 

n(t,r)  = w ( t , r ) + nc(t,r)  (10) 

This  leads  to  the  following  general  representation  of  the 
noise  covariance: 

Kjj  ( t , t 1 , r , r ' ) =Nq6  ( t-t ' ) 6 ( r-r ' ) + Kc  ( t , t 1 , r,r ' ) (11) 

where  Nq  is  the  one-sided  spectral  density  of  the  white 
component  with  units  of  watt/m2-sec.  The  use  of  the  white 
noise  component  to  represent  the  thermal  noise  of  the  array 
is  physically  pleasing  since  the  thermal  noise  is  a broad- 
band noise  which  is  typically  independent  and  identically 
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distributed  from  sensor  to  sensor  in  the  array.  Additionally 
K^(t,t',  r,r')  and  Kp(t ,t ' ,r,r ' ) are  the  covariance  functions 
of  the  total  noise  and  colored  component  respectively. 

The  task  now  is  to  form  both  the  total  noise  and  colored 
covariances.  After  forming  these  covariances,  an  interpre- 
tation of  the  various  jamming  scenarios  and  the  types  of 
covariances  they  imply  is  necessary. 

First,  as  mentioned  before,  a single  jammer  is  represented 
as  a plane  wave  emanating  from  a point  source  in  the  far  field. 
The  complex  envelope  representation  of  the  plane  wave  jammer 
is 


J(t,x,y)  = A(t)exp[j^(t)]exp[j(vxx+vyy)]  (12) 

, COS0  , cos<t> 

where  vx  = — — and  vy  = — 

are  the  spatial  frequencies  of  the  plane  wave  and  0 and  <J) 
are  the  incident  angles  measured  relative  to  the  array's  x 
and  y axes  respectively  (See  Fig  3).  Eq  (12)  can  be  rewritten 
as 


J(t,x,y)  = B(t)exp[j  (v^+Vyy ) ] (13) 

where  B(t)  = A(t)exp[ji|)(t)  ] 

and  B(t)  is  a complex  random  process.  If  more  than  one 
jammer  is  present,  then  the  total  jamming  representation, 
JT(t,x,y),  is: 
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N N 

JT(t,x,y)  = Z J^t.x.y)  = ^B£t)exp[jlvxix+vyiy)]  (14) 

where  N is  the  number  of  jammers  present  and  v . and  v . 
are  the  spatial  frequencies  of  each  individual  jammer.  It  is 
assumed  that  the  jammers  are  uncorrelated,  independent  and 
identically  distributed.  Additionally,  the  expected  value  of 
each  individual  jammer’s  complex  random  amplitude,  B^(t),  is 
zero.  This  leads  to  the  realization  that  the  covariance 
function  of  the  colored  component  is  the  autocorrelation 
function.  The  real  and  imaginary  parts  of  J,p(t,x,y)  have 
the  same  correlation  functions.  It  is  also  assumed  that  the 
jammers  are  temporally  stationary  (recall  Eq  (9a)).  Hence, 
the  general  autocorrelation  function  for  the  N jammers  is: 

Rc(t  ,t*  ,x,y,x’  ,y»  )=E[JT(t,x,y)J*(t',xr,y  ' )] 

N 

= £ Ri(x)exp{j[vlx(x-x' )+viy(y-y' )] } (15) 

where  R.(x)  is  the  individual  temporally  stationary  auto- 
correlation functions  of  the  time  varying  amplitudes  of  the 
jammers.  Note  that  the  individual  terms  in  Eq  (15)  are  in 
the  form  given  by  Eq  (9c)  for  temporally  and  spatially  separ- 
able autocorrelation  functions. 

The  total  noise  autocorrelation  is 

N 

RN(i,r-r' )=No6(x)6(r-r’ )+  E ( x ) exp{ j [ v±x ( x-x ' )+vly (y-y ' ) ] } 
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(16) 


Note  that  the  spatially  colored  component  is  also  spatially 
stationary  due  to  the  properties  of  complex  exponentials. 

This  representation  of  the  total  noise  autocorrelation  is 
very  difficult  to  manipulate  mathematically  in  the  equations 
associated  with  the  signal  detection  problem.  Thus,  to 
further  reduce  this  problem,  consider  representing  R^(x)  as: 

R^x)  = aj_Rj  ( T ) (17) 

This  says  that  all  of  the  jammers  have  the  same  temporal 
statistics  but  may  have  different  average  power  levels.  By 
redefining  the  exponential  in  the  colored  component  as: 

exp[j (vlxx+viyy)]  = $i(r)  (18) 

and 

exp[j  (vixx'+viyy' ) ] = 4>i ( ? ’ ) (19) 

then  Eq  ( 1 6 ) becomes: 

N # 

RN(x,r-r' )=N06(T)6(r-rt  )+Rj(t)  T.  ai<I>i(r)<I>i(r ' ) (20) 

wh  : 

N - * 

Rc(x,r-  / R j ( x ) l ai«J)1(r)<I)1(r ' ) (21) 

Eq  (20)  is  in  the  general  form  given  for  a separable  kernel 
by  Van  Trees  (Van  Trees,  1968:316)  if  the  are  the  eigen- 
values, the  Q.(r)  are  the  eigenfunctions,  and  Rj(t)  is  a 
constant.  If  Rj(x)  is  not  a constant,  the  separable  kernel 
is  still  present  for  only  the  spatial  part  of  the  correlation 
functions.  To  be  able  to  use  this  interpretation,  it  is 
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necessary  to  develop  an  intuition  about  the  kernels  implied 
by  various  jamming  scenarios.  This  intuition  is  based  on  a 
series  representation  of  the  jammers. 

D.  Series  Representation  of  the  J ammlng  Signals 

In  general  the  composite  jamming  signal  can  be  represented 
as  a spatial  Fourier  series  expansion  at  a fixed  time.  This 
expansion  is  given  by: 

00 

J(x,y ) = E J exp[j^Y^x+j^|^  y]  (22) 

n=-« 

where  the  random  variable  Jn  are  given  by 
, L/2  W/2  ~ ~ 

Jn  = LW  1 ’ J(x,y)exp[-j2^!Lx-j— - -y]dxdy  (23) 

-L  -W 
2 2 

When  the  L and  W are  large  enough  the  Jn  are  approximately 
uncorrelated  (Papoulis,  1965:456).  Having  the  uncorrelated 
means  that  exp[j  ^^x+j  ^^-y]  approximates  the  eigenfunctions 
of  a continuous  correlation  function  (Davenport  and  Root, 
1958:97). 

In  general,  a correlation  function  R(r,r')  can  be  repre- 
sented as  a Karhunen-Lo€ve  expansion: 

00 

R(r,r')  = E e.c,(?)c*(?')  (24) 

n=l  1 1 1 

where  the  6^  are  the  eigenvalues  and  the  c^(r)  are  eigen- 
functions of  the  correlation  function  and  a solution  of: 
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1 


L/2  W/2 * 

B.C.(r)  = 1 ! R(r,r* K, (?' )dr'  (25) 

1 x -L  -W  J 

2 2 

An  examination  of  Eq  (2*4)  indicates  useful  interpreta- 
tions implied  by  several  jammer  scenarios.  First,  a single 
jammer  or  a single  jammer  which  is  significantly  more  power- 
ful than  the  other  jammers  approximates  the  condition  when 
there  is  a single  dominant  eigenvalue.  When  there  are 
multiple  jammers  which  are  approximately  equal  strengthed 
and  weak  compared  to  the  thermal  noise  of  the  array  then  all 
of  the  eigenvalues  -re  equal  and  approximately  equal  to  N . 
The  condition  of  multiple  jammers  which  are  of  different 
strengths  and  resolvable  by  the  array  implies  that  the 
eigenvalues  and  eigenfunctions  are  finite  in  number. 

The  last  scenairo  which  was  described  is  the  same 
situation  (finite  number  of  eigenvalues  and  eigenfunctions) 
that  Van  Trees  requires  to  get  his  separable  kernel  f • the 
detection  problem.  Since  it  has  been  shown  that  the  complex 
exponentials  are  the  eigenfunctions  when  (J(x,y)  is  expanded 
such  that  the  Fourier  coefficients  are  uncorrelated,  then 
the  exponentials  are  also  approximately  the  eigenfunctions 
for  the  spatial  correlation  function  of  the  plane  wave 
jammer.  This  spatial  correlation  function  is  given  by  the 
spatial  part  of  Eq  (21): 

N _ * . 

R (r,r')  = Z a.  <J>.  (r)4>.  (r ' ) (26) 

i=l  11  1 
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To  insure  that  the  <t^(r)are  orthogonal  the  following 
criterion  must  be  satisfied 

= sine  (cosG.^  - cosG^  ) ] sine  (cos<J>^-cos<f>j  ) ] (27) 

where  0^  and  0^  are  the  angles  associated  with  x and  x' 
respectively,  and  similarly  for  and  (Jk  and  y and  y'. 

The  term  6. . is  the  Kronecker  delta  function.  The  sine 

ij 

function  is  obtained  by  rewriting  the  exponentials  in  Eq 
(26)  using  Euler’s  equation  and  observing  the  following 
definition  of  the  sine  function 

sine  ( x ) = — --x—  (28) 

Satisfying  Eq  (27)  yields  a finite  set  of  orthogonal  functions 

1 

which  must  be  normalized  by to  get  the  orthonormal  eigen- 
functions. (The  set  of  functions  obtained  visually  from  the 
correlation  function  given  by  Eq  (26)  may  not  be  a complete 
set  of  eigenfunctions  in  which  case  the  set  of  functions 
from  Eq  (26)  would  be  an  approximation  of  the  set  of  eigen- 
functions.) Thus  when  the  <J>^(r)  are  normalized  the  resulting 
correlation  function  is: 

N * 

R ( r-r ' ) = LW  I a, T, ( r ) T, ( r ’ ) (29) 

s i=1  l l l 

4>,(r) 

where  T,  (?)  = and  a,  are  the  eigenvalues.  Eq  (29)  is 

1 /LW  1 

c 
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now  the  proper  form  for  Van  Trees 's  separable  kernel.  Thus, 
when  the  form  of  the  correlation  of  the  N plane  waves  is 
given  by  Eq  (26)  and  the  angular  locations  of  the  jammers 
are  restricted  such  that  Eq  (27)  is  true,  then  the  eigen- 
functions can  be  identified  by  inspection  and  used  to 
develop  the  separable  kernel.  When  the  jammers  are  not 
plane  wave  jammers  or  they  are  not  mutually  orthogonal  then 
the  eigenvalues  and  eigenfunctions  must  be  computationally 
found  by  solving  Eq  (25).  The  difficulty  associated  with 
this  computation  motivates  considering  only  plane  wave 
jammers  who  satisfy  the  criterion  of  Eq  (27)  and  are  thus 
nearly  orthogonal  to  each  other. 

From  the  foregoing  discussion,  the  concept  of  repre- 
senting the  jammers  as  a colored  noise  component  of  the 
noise  in  the  signal  detection  problem  has  been  developed. 
Assumptions  which  facilitate  the  use  of  this  concept  have 
been  established  in  preparation  for  addressing  the  signal 
detection  problem  and  the  performance  evaluation  of  the 
optimum  array  receiver  in  the  presence  of  jamming. 


( 
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Ill . Signal  Detection  Theory 


This  chapter  briefly  reviews  the  formalism  associated  with 
a binary  hypothesis  detection  problem  and  then  identified  the 
general  structure  for  both  the  optimum  array  processor  and 
the  conventional  array  processor.  The  signal  detection  review 
primarily  summarizes  the  results  found  in  Van  Trees  (Van  Trees, 
1968:287-333)  for  a signal  with  single  temporal  and  two  spatial 
dimensions  in  the  presence  of  spatially  non-white  gaussian  noise . 
A.  Review  o£  Signal  Detection  Theory 

The  binary  hypothesis  problem  requires  the  processor  to 
decide  which  of  two  hypothesis  has  occurred  based  upon  the 
received  wave  form.  The  problem  is  expressed  as: 

H-^rCt.x.y)  = s ( t , x ,y ) + n(t,x,y)  te[0,T]  .yet^j^] 

HQ:r(t,x,y)  = n(t,x,y)  (30) 

where  r(t,x,y)  is  the  complex  received  signal,  s(t,x,y)  is  the 
known  signal,  and  n(t,x,y)  is  the  spatially  non-white  gaussian 
noise.  The  noise  term  consists  of  white  and  colored  components 
to  represent  the  presence  of  thermal  noise  and  jamming  signals 
as  previously  discussed.  The  signal  energy  is  given  by: 

T L/2  W/2  * 

E = o -L  -W  s(t>x,y)s  ( t , x ,y )dtdxdy  (31) 

2 2 

where  the  array  is  L by  W and  the  observation  interval  is  [0,T]. 

The  optimum  processor  computes  the  likelihood  ratio  or,  if 
the  likelihood  ratio  is  reducible,  the  sufficient  statistic.  For 
the  purpose  at  hand,  the  sufficient  statistic,  £(r),  is  given  by: 
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(32) 


T L/2  W/2  # 

£(r)  = o _l  r(t,x,y)g  (t,x,y)dtdxdy 
2 2 

where  g(t,x,y)  is  called  the  Fredholm  resolvent  and  is  the  solu- 
tion of  the  integral  equation: 

T L/2  W/2 

s(t,r)  = o l ^ g(t ' ,r* )R  (x,r-r' )dt 'dr ' (33) 

2 -2 

where  R^(x,r-r')  is  the  autocorrelation  function  of  the  noise. 
Fig.  4 illustrates  the  processor. 

Of  course,  once  the  optimum  processor  is  developed,  the 
issue  of  interest  is  measuring  the  performance  of  the 
optimum  processor.  This  measure  is  obtained  by  calculating 
the  probability  of  a decision  error  for  the  processor.  The 
probability  of  error  for  a binary  hypothesis  problem,  where 
the  hypotheses  are  equally  likely,  is: 

P(error)  = erfc(^)  (3*0 

and 

erfc(x)  = -i—  exp  [:i^-]dx  (35) 

x m 2 

where  d is  the  geometric  distance  between  the  means  of  the 
densities  scaled  by  the  standard  deviation  and  given  by  (Van 
Trees,  1968:99): 

j2  ^ [E(£|H1)-E(£|Ho)]2 

d = VarU|HQ)  (36) 

Since  the  erfc(x)  is  a monotonic  function,  a comparison  of 
different  receivers  need  only  involve  the  parameter  d2  (Van 
Trees,  1968:37-38).  The  parameter  d2  is  interpreted  as  the 
output  signal  to  noise  ratio  of  the  optimum  detector  (Adams, 
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1973:16).  This  measure  is  valuable  for  evaluating  receiver 
performance  when  the  signals  are  deterministic. 

B . Optimum  Array 

In  general  to  determine  the  sufficient  statistic  for 
the  optimum  array  detector  it  is  necessary  to  solve  Eq  (33) 
and  then  substitute  the  resulting  Fredholm  resolvent  in  Eq 
(32).  However,  for  a continuous  wave  CW  signal  and  CW 
jammers  there  is  no  time  dependence  in  either  of  these 
signals.  In  other  words,  the  desired  signal  has  constant 
amplitude  as  does  each  of  the  jammers.  Constant  amplitudes 
for  the  jammers  implies  that  the  jammers  temporal  correla- 
tion function  Rj(x)  is  a constant,  Jq . Using  the  CW  signals 
in  Eq  (33)  yields: 

T L/2  W/2 

s(r)  = NQg(t,r)  + JQ  q l W g(t )RS(?-?’ )dt ’dr’  (37) 

-2  -2 

where  R.T(T,r-r')  was  replaced  by  [N  6 ( x ) 6 (r-r ' )+J  R (?-?')]• 

N ’ o os 

Since  there  is  no  time  dependence  in  R (r-r*),  the  time 

S 

integration  is  performed  first  to  give: 

L/2  W/2_ 

s(r)  = NQg(t,r)  + JQ  w g(r'  )Rs(r-r'  )dr ' (38) 

-2  -2 

where 

T 

g(r')  = q g(t ' ,r ' )dt  * (39) 
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There  is  no  time  dependence  in  the  left  hand  side  of  Eq  (38) 
thus  there  can  be  no  time  dependence  on  the  right  hand  side 
of  Eq  (38).  The  only  way  to  realize  the  equality  is  for: 

g(t,r)  = g(r)  (40) 

which  then  means 


g(r' ) = Tg(r ' ) (41) 

Thus  the  Fredholm  resolvent  g(t,r)  has  been  shown  to  be  time 
independent  for  CW  signals  and  CW  jammers.  Eq  (37)  can  now 
be  rewritten  using  this  fact  to  give: 

L/2  W/2 

s ( r ) = NQg(r)  + TJq  g(r' )Rg(r-r ' )dr ' (42) 

2 2 

The  solution  of  Eq  (4l)  for  g(r)  when  substituted  into  Eq  (32) 

yields  the  processor  shown  in  Fig.  5. 

In  general  the  kernel,  H (r-r'),  has  an  infinite  number 

s 

of  eigenvalues  and  eigenfunctions  and  yields  the  following 
solution  for  g(r)  (Van  Trees,  1968:316): 

g(F)  = TJ  E ^5-  C(r)  (43) 

0 i-1  Bi+No 
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where  the  g^  are  the  eigenvalues  and  the  s(r)  are  the  eigen- 
functions of  Rs(r-r').  The  s^  are  given  by: 

L/2  W/2  _ # 

s(r)c1(r)dr  (44) 

2 2 


The  use  of  an  infinite  number  of  eigenvalues  does  not  represent 
a viable  scenario  for  a jamming  environment.  To  get  a viable 
situation  consider  that  R (r-r')  is  a separable  kernel  and  as 
such  has  the  following  solution  for  g(r)  (Van  Trees,  1968:323): 

1 T N Ml  1 

U(F)-TJ  I (45) 

o L i=l  i o J 


For  the  case  of  N nearly  orthogonal  jammers,  it  has  been  demon- 
strated previously  in  Chapter  II  that  for  plane  wave  jammers 
the  spatial  eigenfunctions  and  spatial  eigenvalues  can  be 
obtained  by  inspection  from  the  spatial  autocorrelation  of 
the  jammers.  Thus  using  these  eigenvalues  and  eigenfunctions 
in  Eq  (45)  the  spatial  Fredholm  resolvent  for  CW  plane  wave 
jammers  is: 


g(?)  = i f»(r)-TJ  Ajj-V?)]  (46) 

o l i-1  i o J 

where  T(r)  is  as  defined  in  Eq  (29).  Substituting  this 
definition  into  Eq  (45)  yields: 


g(r) 


s(r)-JQT 


N 

Z 

i-1 


Slal 

VNo 


4'1(r) 

/LW 


(47) 
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where  the  ^(r)  are  as  defined  in  Eq  ( 1 8 ) . 

Once  the  Fredholm  resolvent  is  known,  the  sufficient 
statistic  and  the  performance  can  be  determined  for  the 
optimum  processor.  In  general,  Eq  (36)  can  be  reduced  to 
a simpler  form  for  the  binary  hypothesis  problem  since  the 
noise  is  zero  mean.  The  substitution  of  Eq  (32)  in  Eq  (36) 
yields : 

2 

[T  L/2  W/2  _ * _ 

‘ L W (t,r)dtdr  (lj8) 

-2  -2 J 

uo  ~F  T L/2  W/2  " " " T L/2  W/2 

E o -L  W (t,r)dtdr  £ r*  ( t',  r*)g(  t»,r ' )dt ' dr ' 

2-2  2 2 


which  becomes 


[T  L/,2  Y2  set, 

L°  _?  _p 


r)g*(t,r)dtdr 


d2_  - -2  -? 

o T L/2  W/2  T L/2  W/2 


i 

] 


/ / / 

o L W 

-2  -2 


/ / 
o L 
2 2 


^ H^(x  ,r-r ' )gf  (t  ,r)g(t’,?)dtdt 'drdr ' 


C49) 


Using  Eq  (?"  in  Eq  (*»9)  yields 


T L/2  W/2 


* 


do  = o -L  W s(t,r)g"(t,r)dtdr 
2 -2 


(50) 


which  is  the  general  result  for  the  performance  measure  of 
the  optimum  detector.  For  the  CW  plane  wave  case,  where  the 
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Fredholm  resolvent  is  given  by  Eq  (^5)  and  Eq  (31)  is  true 
for  the  CW  signal,  then  the  performance  measure  is: 


d2 

o 


J T 
o 

N_ 


N 

E 

1=1 


Si2f5i 

No+f3i 


(51) 


The  next  step  is  to  find  the  performance  measure  of  the 
conventional  detector  when  it  is  operating  in  the  same 
jamming  environment. 

C . Conventional  Array  Detector 

The  conventional  array  detector  is  one  which  is  optimized 
for  performance  in  the  presence  of  white  noise  only  (Van 
Trees,  1971:152).  The  development  is  along  the  same  line  of 
thought  as  for  the  optimum  array  detector  except  the  noise 
does  not  have  a colored  component  (which  is  likened  to  saying 
that  the  amplitude  of  the  CW  jammer,  JQ,  is  zero).  With 
this  in  mind,  Eq  (37)  becomes: 


s(r)  = NQg(r)  • (52) 

Thus  the  Fredholm  resolvent  for  the  conventional  array  detector 
is  quite  obviously 

g(r)  = jj-  s(r)  (53) 

o 

The  representation  of  the  conventional  detector  is  commonly 
called  a correlation  receiver.  The  correlation  receiver  is 
also  termed  a beamformer  in  the  case  of  an  array  detector 
(Gallop,  1971:57). 
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Since  the  beamformer  is  optimized  for  white  noise  only, 
it  is  interesting  to  note  its  performance  in  the  presence 
of  the  spatially  correlated  noise  used  to  optimize  the 
optimum  detector . This  noise  is  given  by: 


R^(r-r')  = NQ6(r-r)  + JoRs(r-r') 


(54) 


Substitution  of  Eq  (54)  into  Eq  (49)  yields 


d?  = 


[ 


r L/2  W/2  _ * 

q l y s(r)g  (r)drdt 

-2  -2 


T T L/2  W/2 


T L/2  W/2  T L/2  W/2 


' L J _w  Nng(r)g»  (r)dt'dtdr+  ' / i L _W 


o o -L  -W 
2 2 


o L -W 
-2  2 


-2 


g(r)g  (r ' ) J R (r-r ' ) dtdrdt ' dr ' 
o s 


(55) 


Using  Eq  (52)  in  Eq  (55)  results 


"T  L/2  W/2  # 

* ! 1 s(r)g  (r)drdt 

-L  W 

d*  = 2_-2 

c T T L/2  W/2 

o o -L  -W  s(r)g* (r)drdtdt ' + 
2 2 


in  the  following  expression: 


(56) 


T T L/2  L/2  W/2  W/2 
I I i l I J 

o o -L  -L  -W  -W 

2 2 2 2 


g(?)g  (r')J„R  (r-r ' )dtdt 'drdr ' 

u s 

Substituting  for  R (r-r')  (Eq  (24))  and  having  the  signal 

s 

again  satisfy  Eq  (31),  as  done  for  the  optimum  detector,  then 
the  following  performance  measure  for  the  conventional  detector 
is  obtained: 
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(57) 


! 


Having  the  expression  in  Eq  (57)  permits  comparison  of  the 
conventional  array  detector  with  the  optimum  array  detector 
for  various  jamming  scenarios. 

This  chapter  has  developed  the  tools  necessary  to  evaluate 
the  performance  of  both  the  optimum  and  conventional  array 
detectors.  The  next  chapter  uses  these  tools  to  develop  an 
insight  into  the  performance  improvement  afforded  by  the 
optimum  array  processor  over  the  conventional  processor. 
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IV.  Performance  Improvement 


This  chapter  will  discuss  the  performance  improvement 
that  the  optimum  array  detector  provides  over  the  conven- 
tional beamformer  in  a jamming  environment  composed  of  CW 
plane  wave  signals.  First,  the  method  of  comparing  the 
two  detectors  will  be  discussed.  Next  three  jamming  sce- 
narios will  cover  two  cases  of  multiple  jammers  and  a 
degenerate  case  of  a single  jammer  requiring  only  one 
spatial  dimension  for  processing.  The  two  multiple  jammer 
scenarios  are  first  the  case  of  N nearly  orthogonal  jammers 
and  second  the  case  of  N jammers  which  cannot  be  individu- 
ally resolved  by  the  array.  The  performance  comparisons 
made  for  each  scenario  will  demonstrate  the  superior 
performance  of  the  optimum  array  detector  over  the  conven- 
tional beamformer. 

A . Performance  Measure 

In  order  to  compare  the  performances  of  the  optimum  and 
conventional  detectors  it  is  necessary  to  develop  a measure 
of  the  Improvement.  One  such  measure  is  to  simply  compare 
the  output  signal  to  noise  ratios  for  each  of  the  detectors 
This  comparison  is  denoted  ^ and  is  defined  as 


where  d*  is  the  output  signal  to  noise  ratio  (performance 


measure)  of  the  optimum  array  detector  and  d2  is  the  output 
< signal  to  noise  ratio  of  the  conventional  detector.  Clearly 

the  regimes  of  p which  are  of  greatest  interest  when  p is 
much  greater  than  one  and  when  p is  approximately  equal  to 
one.  The  first  regime  (p  much  greater  than  one)  indicates 
that  the  optimum  array  detector  performs  much  better  than 
the  conventional  detector.  The  second  regime  (p  approxi- 
mately equal  to  one)  indicates  that  both  the  optimum  and 
conventional  detectors  perform  about  the  same. 

Recalling  Eqs  (51)  and  (57),  which  are  tne  general 
expressions  for  the  performance  measures  of  the  optimum 
and  conventional  detectors  in  the  presence  of  CW  plane  wave 
jammers,  the  ratio  p has  the  form: 


E JoT  ? s^i 
N - N~  1 


‘o  1=1  No+gj 


N 

EN  + J T E B.s2 
o o 1 1 


which  can  be  rewritten 


J N 8,  s'  2 

J N 

o *.11 

1 + iT  ,E/iSiZ 

N , .N  +3. 

o i=l  o ^i _ 

o i=l 

where 


( 


* 


£7?  W72 
\ \ 
-L  -W 
2 1 


_ l s(x,y)s  (x,y)dxdy 


(59) 
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This  normalizes  the  signal  energy  and  leads  to  an  expression 
which  is  more  easily  interpreted.  The  expression  in  Eq  ( 6 0 ) 
is  valid  for  comparing  the  performance  of  the  optimum  and 
conventional  detectors  in  an  environment  of  CW  jammers  so  long 
as  the  colored  noise  component  results  in  a separable  kernel. 
The  discussions  are  all  special  cases  of  Eq  (60). 

B.  N Nearly  Orthogonal  Jammers 

The  scenario  involving  N nearly  orthogonal  CW  plane  wave 
jammers  requires  that  the  criterion  established  in  Eq  (27)  be 
met.  With  this  criterion  met,  Eq  ( 6 0 ) becomes: 


p = 1 + (LW)2J: 
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where  the  3^^  are  LWa.^  (recall  from  Eq  (29)  and  the  s^  are 
given  by: 


..  L 

w - 

sin(vsx  vnxl)2 

sin( vsy-vny i ) 2 

_l(vsx“vnxi) 

_ ^jy-Vi5 

(64) 


( 


33 


where  vgx,  vgy,  vnxl,  and  vnyi  are  the  spatial  frequencies 
of  the  signal  and  jammer  and  are  given  by: 
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nyi  A 


(65a) 


( 65b ) 


( 65c ) 


( 6 5d ) 


n 


and  the  0s,  0ni,  <j>s  and  <J>nl  are  the  incident  angles  rela- 
tive to  the  array's  x and  y coordinate  (recall  Fig.  3)  for 
the  signal  and  jammers.  The  eigenfunctions  used  to  compute 
the  s^  are  given  by  Eq  (29)  as  previously  discussed.  After 
applying  the  assumption  that  the  signal  and  jammers  are  at 

the  same  wavelength,  A , and  recalling  the  definition  of 

s 

sinc(x)  in  Eq  (28),  Eq  (64)  becomes 


3'^  = sine  £^(cos0g-cos0n^)  J sinc£j^(cos<(>s-cos<|)nl) J (66) 


Thus  the  significant  aspect  of  Eq  (63)  with  regard  to  the 
performance  improvement  is  the  value  of  the  s^2. 
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The  double 


sum  coefficient  of  the  s’^2  in  Eq  (63)  is  always  a positive 
number  and  serves  to  weight  the  terms  involving  the  s'^2. 
Several  observations  can  be  made  relative  to  the  performance 
improvement . 

First,  when  all  of  the  s'^2  are  zero,  then  the  performance 
indicator,  y is  unity.  This  situation  occurs  when  the  signal 
is  nearly  orthogonal  to  all  of  the  jammers.  In  this  case, 
the  conventional  detector  can  detect  the  signal  as  well  as 
the  optimum  detector  since  the  effect  of  the  jammers  has  been 
neutralized  by  the  mutual  orthogonality  of  the  signal  and 
jammers. 

Second,  y is  also  unity  if  the  s'^2  all  sum  to  one.  This 
situation  occurs  if  the  signal  is  perfectly  aligned  with  one 
of  the  jammers.  Then,  in  this  case,  both  the  optimum  and 
conventional  detectors  perform  equally  poorly. 

The  above  discussion  addressed  the  cases  of  equal  per- 
formance for  the  optimum  and  conventional  detectors.  These 
cases  are  interesting,  but  the  case  of  greatest  interest  is 
that  in  which  the  optimum  detector  is  clearly  superior  to  the 
conventional  detector  (y  is  much  greater  than  unity).  To 
look  at  this  case,  consider  first,  that  due  to  the  mutual 
orthogonality,  the  signal  will  project  only  on  those  jammers 
who  fall  within  the  resolution  range  of  the  array  when  it's 
main  beam  is  directed  at  the  signal.  All  other  jammers  will 
have  no  effect  upon  the  signal. 


For  the  case  when  there  are  two  jammers  close  enough 
for  the  signal  to  project  on  each  of  them,  the  Eq  (63) 
becomes : 
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Eq  (67)  can  be  simplified  greatly  if  the  following  realtion- 
ships  are  assumed: 


“l  ” a2 


(68a) 


-12  _ Ql  2 

S1  ' S2 


(68b) 


Eq  (68a) implies  the  two  jammers  are  of  equal  strength.  Eq 
(68b)  implies  that  the  signal  projects  equally  onto  each 
Jammer.  Stated  another  way,  the  signal  is  angularly  equi- 
distant from  each  jammer  in  either  the  x or  y dimension. 


Then  Eq  (67)  becomes 
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or  by  letting  y be  given  by 
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Y 


(70) 
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then  Eq  (69)  becomes 


y = 1 + Y [s,12(|  - s^2)]  (71) 

where  s^2  is  given  by 

s^2  = sinc2^^  (cos0 

From  Eq  (72)  and  the  assumptions  made  in  Eq  (68b),  it  is 
clear  that  s!^2  is  bounded  by  zero  and  one-half.  Fig.  6 plots 
Eq  (71)  as  a function  of  s^2.  From  Fig.  6,  the  maximum  value 
of  y occurs  when  s^equals  one-fourth  and  is  given  by 

>w  - 1 + -°62^  (73) 

This  situation  represents  the  fact  that  the  signal  projects 
twenty  five  percent  of  its  power  in  the  direction  of  each 
of  the  two  jammers.  For  y to  be  much  greater  than  unity, 
it  is  necessary  for  the  following  condition  to  be  true: 


-cos6  . ) 
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(72) 


Ys^2(|-  - s^2)>>  1 


(74) 


PERFORMANCE  MEASURE 


Recalling  Eq  (69),  the  condition  described  in  Eq  (7*0  occurs 
when  the  following  situation  exists: 
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(75) 


which  says  that  twenty  five  percent  of  the  thermal  noise  (Nq) 
to  interference  power  (J^a^LW)  ratio  is  much  less  than  the 
signal  power  in  the  interference  direction.  Thus  by  knowing 

\ 

the  thermal  noise  of  the  array,  the  array  size,  the  inter-  j 

ference  power  and  the  /location  of  the  jammers  in  both  spatial  I 

J 

dimensions,  the  signal  location  can  be  positioned  to  achieve  j 

i 

the  criterion  established  in  Eq  (75).  Alternately,  knowledge 

< 

and  ability  to  change  any  of  the  parameters  permits  the 

( 

optimum  detector  to  be  peaked  for  maximum  performance  improve- 
ment for  a given  situation. 

As  an  example,  consider  letting  the  ratio  of  the  inter- 
ference power  to  the  thermal  noise  be  30dB  and  the  value  for 
s^2  be  0.1,  then  the  criterion  in  Eq  (75)  are  satisfied  and 
substitution  into  Eq  (69)  Indicates  that  the  performance 
improvement  of  the  optimum  detector  is  on  the  order  of  68dB. 

Computing  d2  using  Eq  (51)  and  d2  using  Eq  (57)  for  the  same 
condition  yields: 

d2  = lOldB  (76a) 

d2  = 33dB  (76b) 
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These  results  are  for  a normalized  signal  energy  and  indicate 
that  the  optimum  detector  is  clearly  better  than  the  beam- 
former  . 

Before  concluding  the  mutually  orthogonal  jammer 
scenario,  an  interesting  degenerate  case  is  the  one  for  only 
a single  jammer  present.  For  the  single  jammer,  Eq  (63) 
becomes : 


I"  ( LW ) 2 a 2 J 2 ‘ 
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where  s'  2 is  given  by  Eq  (72)  when  the  subscripts  are  deleted. 
For  the  single  jammer  case,  s'  2 is  bounded  by  zero  and  one. 
Fig.  7 plots  the  following  function  versus  s'  2: 

M = 1 + y'  s'  2 (1-s'  2 ) (78) 

where  y'  is  the  constant  coefficient  in  Eq  (77). 

From  Fig.  7,  the  maximum  value  occurs  when  s'  2 equals 
one-half  and  is  given  by: 

“max  - 1 + ' 2 5V  (79) 

This  situation  represents  the  fact  that  the  signal  projects 
fifty  percent  of  its  power  in  the  direction  of  the  jammer 
(Pasupathy,  1978:161).  For  y to  be  much  greater  than  unity 
it  is  necessary  for  the  following  criterion  to  be  true: 
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(80) 


This  condition  occurs  when 


J oLW 
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which  says  the  ratio  of  thermal  noise  to  interference  power 
is  much  less  than  the  signal  power  in  the  interference  direction. 
Using  the  same  example  as  used  previously  for  the  two  nearly 
orthogonal  jammers,  the  performance  improvement  of  the  optimum 
detector  over  the  beamformer  is  now  60dB.  Again  the  example 
was  done  for  a unit  energy  signal. 


PERFORMANCE  MEASURE 


C.  N Non-Resolvable  Jammers 


The  next  jamming  scenario  of  interest  is  the  scenario 
in  which  the  array  cannot  resolve  the  individual  jammers 
due  to  insufficient  angular  separation  among  the  jammers. 
This  case  corresponds  to  having  a single  spatially  band- 
limited  jammer  with  N significant  eigenvalues.  If  the  N 
jammers  are  all  of  equal  power  then  the  eigenvalues  can  be 
considered  to  be  approximately  equal.  Thus  the  expression 
for  p given  by  Eq  ( 6 0 ) becomes: 
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where  n is  the  significant  eigenvalue  and  the  s'^  are  as 
previously  defined  in  Eq  ( 6 1 ) and  Eq  (4*0  . The  notable 
difference  between  this  case  and  the  case  of  N nearly 
orthogonal  jammers  is  the  determination  of  the  eigenfunctions. 

For  the  orthogonal  jammers,  the  eigenfunctions  were  the  weighted 
exponential  propagation  functions  of  each  jammer.  For  the  band- 
limited  case,  the  eigenfunctions  need  to  be  computed.  This 
determination  of  the  eigenfunctions  must  be  done  for  specific 
examples  and  cannot  be  done  in  general.  However,  by  inter- 
preting the  summation  of  the  s^2  as  the  signal  power  in  the 
interference  beamwidth  of  the  array  (Pasupathy,  1978:161), 
it  is  possible  to  develop  an  intuition  about  the  performance 
improvement  of  the  optimum  detector  versus  the  conventional 
detector.  By  redefining  the  summation  of  the  s'^2  Eq  (82)  becomes 
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Comparing  Eq  (83)  with  Eq  (77),  it  is  obvious  that  the  two 
equations  are  of  the  same  form,  differing  only  in  the  constant 
coefficient.  Thus,  the  same  analysis  is  valid  as  was  performed 
for  the  single  jammer  in  the  previous  section  with  some  ampli- 
fying observations  regarding  the  physical  interpretations. 

When  o2  is  equal  to  one-half,  p achieves  its  maximum  value 
and  is  given  by: 


when  n is  much  greater  than  the  thermal  noise,  NQ.  To  inter- 
pret Eq  (85),  it  is  necessary  to  look  at  d2  and  d2 . These 
quantities  are  given  by 
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Evidently,  d2  is  independent  of  n for  the  given  situation 
but  d2  depends  inversely  on  the  relationship  between  pJQ 
and  N . Thus  the  performance  improvement  indicated  by 
Pmax  primarily  due  to  a degradation  of  the  conventional 
detector  performance  rather  than  an  increase  in  the  perfor- 
mance of  the  optimum  detector  (Pasupathy,  1978:162).  This 
was  not  the  case  for  the  orthogonal  jammer  scenario.  In 
that  scenario  the  performance  improvement  was  tied  to  both 
the  improvement  afforded  by  the  optimum  detector  and  the 
degradation  of  the  conventional  detector. 

To  quantify  this  scenario  in  terms  of  performance  improve- 
ment as  previously  done,  consider  letting  r|J  be  much  greater 

pj  ° 

than  Nq,  and  the  ratio  - be  equal  to  30dB,  then  for  a a2 

o 

equal  to  0.1  the  performance  improvement  is  approximately 
19.5dB  when  all  of  the  assumptions  and  conditions  are 
substituted  into  Eq  (83).  The  fact  that  this  scenario  has  N 
jammers  is  contained  in  the  a2  and  variations  in  the  results  due 
to  a different  number  of  jammers  are  attributed  to  variations  in 
a2. 

D.  Single  Jammer,  Single  Spatial  Dimension 

The  final  jamming  scenario  is  really  a degenerate  case  of 
the  single  jammer  scenario  discussed  for  two  spatial  dimen- 
sions. This  degenerate  case  is  one  in  which  it  is  necessary 
to  look  only  at  a single  spatial  dimension.  The  degenerate 
case  comes  about  when  either  of  the  sine  functions  in  Eq  (66) 
is  unity.  This  means  that  the  signal  and  jammer  are  aligned  in 
one  of  the  spatial  dimensions.  With  this  being  the  case,  it 


Y 


is  necessary  to  look  only  at  the  single  spatial  dimension 
to  determine  the  performance  improvement.  Thus  for  the 
one  dimension  case  Eq  (77)  becomes: 
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where  q2  is  given  by: 


q2  = sine 2 [^(  cos0  -cos©  )] 
a s n 
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(88) 


Again  the  same  interpretations  as  made  for  the  two  dimension 
single  jammer  case  apply.  Fig.  8 illustrates  a specific 
example  of  the  one  dimension  case.  Given  that  a,  Nq,  As,  and 
L are  all  fixed  for  the  example,  it  is  easy  to  see  the  effect 
of  different  angular  separations  between  the  signal  and  the 
jammer.  Thus,  Fig.  8 indicates  the  improvement  for  different 
signal  incident  angles  (0  ) to  the  array  of  the.  optimum 

O 

detector  over  the  conventional  detector.  For  example,  given 
that  the  array  could  be  easily  reoriented  so  as  to  change  the 
incident  angle  of  the  signal,  then  knowledge  of  the  angular 
difference  (A0)  between  the  signal  and  jammer  would  indicate 
which  incident  angle  would  provide  the  greatest  performance 
Improvement.  Also  from  Fig.  8 the  performance  improvement 
for  the  example  which  has  been  used  in  all  of  the  previous 
scenarios  is  given  to  be  approximately  20dB. 

This  chapter  has  succeeded  in  developing  a performance 
measure,  p , which  was  used  to  evaluate  three  jamming  scenarios. 
The  three  jamming  scenarios  were  evaluated  with  a common 

c 

45 


it  5 
if.l 


Figure  8.  Angular  differ 


example  and  It  was  shown  that  the  optimum  array  detector 
is  significantly  better  than  the  beamformer.  Additionally 
physical  interpretations  were  made  to  enhance  the  informa- 
tion provided  by  the  performance  comparisons.  The  final 
conclusions  of  this  thesis  and  recommendations  for  additional 
research  are  addressed  in  the  next  chapter. 
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V.  Conclusions  and  Recommendations 


A . Conclusions 

The  preceeding  discussion  has  developed  three  areas. 

First,  the  representation  of  several  jamming  scenairos  as 
a colored  noise  component  in  the  noise  term  of  the  binary 
hypothesis  detection  problem.  Second,  a measure  of  perfor- 
mance improvement,  y,  was  developed  by  comparing  the  output 
signal  to  noise  ratios  of  the  optimum  and  conventional  array 
detectors.  Analysis  of  the  performance  improvement  was 
presented  for  the  jamming  scenarios  of  interest  based  on 
total  knowledge  of  the  desired  signal  and  the  jamming  signals. 
This  analysis  developed  relationships  which  identified  the 
parameters  which  the  system  designer  could  address  to  achieve 
maximum  improvement  from  the  optimum  array  detector.  Finally, 
an  example  was  presented  for  each  scenairo  which  showed  that 
the  optimum  array  detector  was  significantly  better  than  the 
beamformer . 

Several  of  the  points  previously  made  need  to  be  reiterated. 
First,  for  single  Jammers,  the  maximum  performance  improvement 
is  achieved  by  the  optimum  detector  when  fifty  percent  of  the 
signal  power  is  projected  in  the  jammer's  direction.  This  is 
also  valid  for  multiple  jammer  situations  when  the  signal  is 
nearly  orthogonal  to  all  but  one  of  the  jammers.  Secondly, 
for  single  Jammers  and  N nearly  orthogonal  jammers  the  thermal 
noise/interference  power  ratio  must  be  much  less  than  the 
signal  power  in  the  jamming  direction  to  achieve  significant 
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performance  improvements.  In  the  case  of  the  non-resolvable 
jammers  it  was  seen  that  significant  performance  improvement 
was  largely  due  to  degradation  of  the  conventional  detector 
rather  than  the  improved  performance  of  the  optimum  detector. 

These  observations  bring  out  two  points  which  are  appli- 
cable to  the  array  problem.  First,  the  signal  power  and  the 
projection  of  the  signal  power  in  the  interference  (jamming) 
direction  significantly  affect  the  performance  of  the  detectors. 
Secondly,  the  thermal  noise,  in  the  case  of  the  non-resolvable 
jammers  particularly,  is  a significantly  contributing  factor 
in  the  performance  of  the  detectors. 

It  is  recognized  that  the  criterion  which  allows  the 
single  spatial  dimension  degeneration  is  probabilistically 
very  unrealizable.  None  the  less,  the  single  spatial  dimension 
provides  the  ability  to  get  a "back  of  the  envelope"  approxi- 
mation which,  when  compared  to  the  two  spatial  dimension 
results,  provided  the  engineer  with  an  insight  into  the  problem 
without  a great  deal  of  effort. 

Finally,  a single  example  was  exercised  for  each  secnario. 

In  this  example  the  ratio  of  jammer  interference  power  to 
thermal  noise  was  equal  to  30dB  and  the  parameter  related  to 
the  projection  of  the  signal  onto  the  jammer  was  equal  to  0.1. 

The  resulting  performance  improvements  were  68dB,  60dB,  19.5dB, 
and  20dB  for  the  two  orthogonal  jammers,  single  jammer  (two 
spatial  dimensions),  spatially  bandlimited  jammer,  and  the  single 
Jammer,  (single  spatial  dimension) respectively . The  variations 
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in  the  performance  improvement  are  related  to  the  complexity 
of  the  individual  jamming  scenario.  However,  all  of  the 
scenarios  indicate  that  the  optimum  array  detector  is  clearly 
better  than  the  beamfcrmer. 

B.  Recommendations 

The  previous  work  is  based  on  full  knowledge  of  the  signal 
and  the  jammers,  particularly  with  regard  to  amplitude  and 
direction.  This  fact  suggests  that  future  research  should 
extend  the  assessment  of  performance  improvement  by  relaxing 
the  assumptions  which  were  made  in  this  study.  Specifically, 
addressing  the  performance  improvement  for  time  varying 
amplitudes  for  both  the  desired  signal  and  jammers  and  the 
inclusion  of  temporally  dependent  signals  would  be  the  next 
logical  extension  of  this  research.  Succeeding  extensions 
of  the  effort  would  include  examination  of  signals  with 
random  amplitude  and  random  phase,  and  a relaxation  of  the 
directional  knowledge  of  the  signal's  location.  Completion 
of  these  research  tasks  would  provide  a relatively  complete 
picture  of  the  performance  improvement  that  an  optimum  array 
detector  can  provide. 
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